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ABSTRACT 


Strain energy release rates for edge -delaminated composite 
laminates were obtained using quasi three-dimensional finite element 
analysis. The problem of edge-delamination ap the -35/90 interfaces of an 
eight ply [0/±35/90] s composite laminate subjected to uniform axial strain 
was studied. A quasi three-dimensional finite element analysis was used to 
calculate the total and individual components of the strain energy release 
rate. The individual components did not show convergence as the 
delamination tip elements were made smaller. In contrast, the total strain 
energy release rate, G, converged and remained unchanged as the delamination 
tip elements were made smaller and agreed with the total G calculated using 
a closed- form equation derived from the rule of mixtures and classical 
laminated plate theory. The studies of the near -field solutions for a 
delamination at an interface between two dissimilar isotropic or orthotropic 
plates showed that the imaginary part of the singularity is the cause of the 
nonconvergent behavior of the individual components . To evaluate the 
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accuracy of the results, an eight ply [0/±35/r/90] s laminate with the 
delamination modeled in a thin resin layer, similar to the resin layer that 
exists between the -35 and 90 plies, was analyzed. Because the delamination 
exists in a homogeneous isotropic material, the oscillatory component of the 
singularity vanishes. The strain energy release rates remained unchanged as 
the delamination tip elements were made smaller. Comparison of the strain 
energy release rates for the 'bare' interface laminate, i.e. one without the 
resin layer, and for the laminate with the resin showed that the 'bare' 
interface models are a very good approximation for the resin case if the 
delamination tip elements were one-quarter to one-half of the ply thickness. 


INTRODUCTION 

Composites are used extensively in aerospace, automobile, and civil 
engineering structures because of their high strength- to -weight ratios. 
Because delamination is a common failure mechanism of composite laminates, 
research has been directed toward understanding delamination mechanisms and 
quantifying the strength of composite laminates with delaminations [1,2]. 

Recently, to quantify the interlaminar fracture toughness of 
composite laminates, the edge-delamination test was proposed in references 
1-3. For the family of [O/±0/9O] s laminates subjected to remote uniform 
axial strain, O'Brien [2] found that for a given interlaminar fracture 
toughness the minimum failure strain was obtained when 9 - 35°. Therefore, 
the [0/±35/90] s family of laminates, consisting of [±35/0/90] s , 

[35/0/-35/90] s , and [0/±35/90] s layups, was recommended for the edge- 
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delamination test. In this family of laminates, the interlaminar normal 
stresses are greatest when the delamination is between the 90° ply and its 
neighbor. For these laminates, the total strain- energy release rate for 
self-similar delamination growth can be obtained using an equation derived 
from the rule of mixtures and the classical laminated plate theory (CLT) 
[1], However, the individual components (mode I, mode II, and mode III) 
cannot be determined from simple formulas like that given in reference 1. 
To determine the individual components, the boundary value problem of the 
laminate with edge -delaminations needs to be solved. The quasi -three 
dimensional (Q3D) finite element analyses are useful to determine the 
interlaminar stresses and the strain energy release rates for delamination 
growth . 


The laminates used in the edge -delamination test are cocured. 
Experimental evidence indicates that very thin resin- rich layers of about 
0.0004 in. thick exist between neighboring plies [5,6]. The Q3D finite 
element analysis should model these resin layers as well as the individual 
plies. Because the resin layers are very thin and smooth transitions are 
needed in the model, large numbers of elements and nodes are required to 
model the resin layers. Therefore, the resin layers are usually neglected 
and the delamination is assumed to be at a discrete interface between 
neighboring plies [1-4,7]. In reference 3, a [0/±35/90] s laminate with a 
'bare' interface was modeled with 8-noded isoparametric parabolic elements 
everywhere. At the delamination tip, square non-singular parabolic elements 
were used. The convergence of the strain energy release rates was studied 
for various size delamination tip elements. The study showed that the 


3 



individual components of the strain energy release rates did not converge as 
the size of the delamination tip elements was reduced. Therefore, the 
accuracy of the individual components determined by such a model is 
questionable. In contrast, the total strain energy release rate converged to 
show no change with the reduction of the size of delamination tip elements . 
Also the total strain energy release rate agreed extremely well with that 
calculated using the CLT formula from reference 1. 

The first objective of this paper is to identify the source of the 
non- convergent behavior of the individual components of strain energy 
release rate. The second objective is to establish the accuracy of the 
finite element solution for the ’bare' interface problem. To this end, a 
cocured laminate with res in- rich interfaces was analyzed and compared to a 
'bare' interface model. 

First, the edge -delamination problem and the Q3D finite element 
analyses are reviewed. Then the strain energy release rates for the 
[0/±35/90] s laminate with a 'bare' interface between -35 and 90 plies were 
obtained using either non-singular or quarter-point singularity elements at 
the delamination tip. The singularity elements were used to determine if 
the use of the non- singular elements was the cause for the non -convergent 
behavior of the individual components. While the quarter-point elements 
produce the classical square root singularity at the delamination tip, the 
singularity at the delamination tip between two dissimilar materials is 
known to be of the form -1/2 ± iY for two-dimensional (2D) problems [7-11]. 
The imaginary part of the singularity, the so called "oscillatory" 
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component, y , may also be a cause for the non- convergence of the individual 
components. To study this aspect, the near field stress distribution for 
cracks between dissimilar isotropic or ortho tropic materials under plane 
stress or plane strain was examined. Next, the edge -delamination problem 
with delamination between two different isotropic materials was studied 
using the Q3D finite element analysis with the oscillatory component either 
zero or non zero. Last, the strain energy release rates for a [0/±35/90] s 
laminate obtained with a resin interface are presented for various size 
delamination tip elements. These results are compared with those obtained 
without modeling the resin layer. 
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SYMBOLS 

length of delamination 
half width of the laminate 

Young's modulus for orthotropic material in the j- 
direction 

Young's modulus of resin 
total strain energy release rate 

mode I, mode II, and mode III strain energy release rate 

components , respectively 

shear modulus for orthotropic material 

ply thickness 

resin thickness 

Pi 

total thickness of the laminate 
displacement functions 


5 


u,v,w 


x,y,z 


A 

G o 

{e} 


{a} 

y 

v j 

V jk 

v r 

0 


displacements in the x- , y- , and z-directions , 

respectively 

Cartesian coordinates 

delamination tip element size 

uniform axial strain in the x- direction 

Cartesian engineering strains, { e x , e y , e z , e xy , e yz , 

G zx 1 

shear modulus for the j** 1 isotropic material 
Cartesian stresses, { a x , a y , o z , <^ X y» a yz » a zx ) 
oscillatory part of the singularity 
Poisson's ratio of the j isotropic material 
Poisson's ratio for orthotropic material 
Poisson's ratio of the resin 

angle between x-axis and the fiber axis (see fig 1(a)) 


Superscript 

bar denotes complex conjugate 


EDGE -DELAMINATION ANALYSIS 

A symmetric eight ply laminate subjected to uniform axial strain e G in 
the x-direction is shown in Fig. 1. Delaminations at both edges are located 
symmetrically about the laminate midplane (at the -35/90 interface above the 
midplane and at the 90/-35 interface below the midplane). Because the 
laminate is long in the x-direction, all x—constant planes away from the 
ends deform in the same manner. Therefore, away from the ends the 
displacements are assumed to be, 
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( 1 ) 


u- U(y,z) + e Q x 
v- V(y,z) 
w- W(y , z) 

where U,V,W are displacement functions expressed in terms of y and z alone 
[4,7,12,13]. Equations (1) describe a 'quasi three-dimensional' (Q3D) 
problem. The modifier 'quasi' is used because there are three displacements 
in the three directions, but the gradients of U, V, and W with respect to 
the x-coordinate are zero. Thus only an x=constant plane needs to be 
analyzed to obtain the stresses in the laminate. Figure 1(b) shows a 
typical x-constant plane. Because of symmetries, only one quarter of this 
plane was analyzed (see Figure 2) . 

The material properties used in this study were 

E n - 19.5 x 10 6 psi ; E 2 2 “ E 33 - 1.48 x 10 6 psi 

v 12 “ v 13 “ ° • 3 I v 23 “ - 49 

G 12 - G 13 - 0.8 x 106 psi ; G 2 3 - 0.497 x 10 ^ psi 

where the subscripts 1, 2, and 3 correspond to the longitudinal, transverse, 
and thickness directions of a zero degree ply [3,7]. The resin properties 
used in this study are E r - .5 x 10^ psi ; v r - . 3 . These properties gives 
a shear modulus of 0.192 x 10^ psi. 

When subjected to axial strain, the 0/±35 ply group above the 
delamination tends to contract (and shear) by different amounts in the y- 
direction (and the xy-plane) compared to the 90 ply below the delamination. 
Therefore, to maintain the displacement continuity along the -35/90 interface 
(z-h in Figure 2) , Oy and a X y stresses develop in the interior as shown in 
figure 2. Equilibrium requires that stresses Oy Z , a xz , and o z exist along 
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the z - h interface (see fig 2(b)). Because the faces of the delamination 
are stress free, these interlaminar stresses exist only in the range 0 < y < 
(b-a). These interlaminar stress give rise to three modes of deformation at 
the delamination tip. 

Figure 3(a) shows modeling for the edge -delamination in Fig 1(b). Two 
types of modeling were used. In the first model, 8-noded isoparametric 
parabolic elements were used to model the region. In this model, these non- 
singular elements were used everywhere including the delamination tip. The 
second model was identical to the first, except collapsed quarter-point 
singularity elements were used around at the delamination tip, as shown in 
Fig.. 3(b). The details of the Q3D analysis are given in reference 4 and 
hence, are not repeated here. 


The total strain energy release rate, G, was obtained from CLT as, 


f e LAM - E ] 


( 2 ) 


2n 

where E lam and E are axial stiffnesses calculated from the classical 
laminated plate theory (CLT) for the undelaminated and completely 
delaminated laminates (along one or more interfaces), respectively [1,2]. 
The n is the number of delaminated interfaces , and t is the total laminate 
thickness. As previously mentioned, the individual components of strain 
energy release rate cannot be obtained by simple formulas like equation 2. 
However, these components can be calculated using a Q3D finite element 
analysis [1-4,7]. 
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The components of strain energy release rate were calculated using 
forces and displacements near the delamination tip with Irwin' s virtual 
crack closure technique (VCCT) as , 

1 

Gi - - [ F zi { t n (W m - W m ') + t 12 (W k - W k ') } 

2A . (3) 

+F zj { t 21 (W m - W m ') + t 22 (W k - W k ') } ] 

where F z ^ is the force in the z-direction at node i, W m is the W- 
displacement at node m, etc. (see Fig. 4) [14,15]. For non-singular 
elements (see Fig 4(a)), t^i - 1 ; t^ 2 ** 0 ; t 2 ^ * 0 ; t 22 " 1 and for 
quarter-point singularity elements (see Fig 4(b)), t^i - 6-(3 tt/ 2) ; t^ 2 — 
6 tt- 20 ; t 2 ^ - 1/2 ; t 22 — 1. Similar expressions for Gjj and Gjjj were used 
with F z replaced by Fy and F x and with W replaced by V and U, respectively. 

Several finite element idealizations were used with both the non- 
singular and singular elements. Delamination tip element sizes were A/h - 
0.5, 0.25, 0.125, and 0.0625. In the finite element analysis with non- 
singular elements, the singularity at the delamination tip is not modeled. 

On the other hand, when the quarter-point elements are used the square root 
singularity is incorporated at the delamination tip. The VCCT was used with 
both element types to evaluate the strain energy release rates. Table 1 
presents the individual and the total strain energy release rates calculated 
with non-singular and quarter-point singularity elements, respectively. 

These results are also plotted in Fig 5. Both the non-singular and the 
singular elements yielded similar results. They showed that the total 
strain energy release rate remained unchanged as the size of the 
delamination tip elements decreases. However, both element type shows a 
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slight non- convergence for the individual modes. This situation is in 
contradiction to the conventional belief that better accuracy may be 
obtained when the size of the crack/delamination tip elements are decreased. 
The mode -I strain energy release rate becomes larger as the delamination tip 
element becomes smaller (i.e. for smaller values of A/h) . The mode -II 
strain energy release rate shows the opposite trend. The mode -III component 
has a very small negative value. The total strain energy release rate 
computed from the finite element analysis agrees very well with that 
calculated by equation (2) . 


ANALYSIS OF NON- CONVERGENCE 

In this section, the reasons for the non- convergence of the 
individual components of the strain energy release rate are explored. In the 
edge -delamination problem the delamination exists at the interface between 
two different anisotropic materials. For this case the singularity is not 
the classical square root singularity but is of the form , where r 

is the radial distance measured from the delamination tip. The 7 depends 
on the material properties of the two materials [7-11]. This imaginary 
power leads to the stress oscillations very close to the delamination tip. 
The oscillatory component of the singularity may cause the non- convergence 
of the individual G components. To explore this possibility, a simpler 
problem of a bi-material plate in plane stress or plane strain was studied. 

The problem of an interfacial crack in a bi -material plate with two 
different isotropic materials in either plane stress or plane strain was 
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examined. Then, the same problem with two orthotropic materials was 
studied. 


Interfacial Crack in Bi-Material Isotropic Plate 

Figure 6 shows a bi -material plate with an inter facial crack 
subjected to remote uniform tension. Due to the applied stress, each of the 
two materials tries to contract differently. Because they are joined at the 
interface (z-0) , to maintain continuity a shear stress <jy Z develops along 
the bond line. Thus a mixed mode condition develops for remote tension. 
However, for certain combinations of material properties that satisfy Eq. A5 
(see appendix) , the contractions very close to the crack tip are identical 
and hence the shear stresses Oy Z are not required to maintain the 
continuity. On the other hand, if two materials do not satisfy Eq. A5, then 
shear stress stresses Oy Z are needed on the z-0 line to maintain 
compatibility and the near field stress state will be different from the 
classical square root distribution and is of the form, 

o z - ( Pr"- 1 -/^ + Pr"*-/^ -i?) + . . . 

Oy Z - ( Rr' 1 / 2 +iy + Rr -1 / 2 ' iY ) + ... (4) 

where P and R are complex constants, and P and R are their corresponding 
complex conjugates [8-11]. These constants depend on the loading. The 
oscillatory power, V , in Eq. (4) is given by, 

1 + Mb k A 

y In [ ] (5) 

2jt kg 

where /ij are the shear moduli and kj are given by 
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for plane strain and 


and j— A,B. 


3 - 4 v. 


3 - 


1 + v i 


for plane stress 


When the two materials satisfy Eq. A5 the oscillatory power, y , is 
identically zero because the terms in the square brackets in Eq. 5 become 
unity. 


From the results in reference 11, it can be shown that the 
constants P and R in Eq. 4 are related to each other as R - -iP . The 
relative displacements of the crack faces behind the crack tip have the form 

v «* ( + Sr'^/^-iy ) + . , . 

( 6 ) 

W - ( QrV2+iV + q r -l/2-iy ) + ... 

where the constants Q and S are complex constants and are related to each 
other as S - -iQ. 

The strain energy release rates for the crack along the interface of 
this bi-material plate can be determined as outlined below. 

Strain Energy Release Rates 

The strain energy release rates were obtained by using the near 
field solution of Eqs. 4 and 6 and Irwin's VCCT as, 
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( 7 ) 


1 A 

G - lim / [ a z (r) .w(A-r) + a yz (r) .v(A-r) ] dr 

A-0 2A 0 


where A is now the distance over which the crack is assumed to close. (Note 
that in the finite element analysis A was the delamination tip element 
size.) 


The mode I and mode II components of the strain energy release rates 

are 




1 

A 



G I 

- lim 

— 

J 

[ a z (r).w(A-r) ] dr 



A-0 

2A 

0 





1 

A 



G II 

- lim 

— 

/ 

[ c7y Z (r) .v(A-r) ] dr 

(8) 


A-0 

2A 

0 



Substituting the stresses and displacements from Eqs 4 and 6 into Eqs. (8) 
gives , 


G! - lim [ PQA 2iY I* + PQA' 2iY Ii + PQ I 2 + PQ I 2 ] 

A-0 (9) 

where 

n/2 

11 - J cos 2 /} (sin/3 cos£) 2 *- y d/3 

0 

n/2 (10) 

12 - J cos 2 /) (cos/9/sin/9) 2 * Y d/3 

0 

- (ir/2) (2iY +!)/( e” y + e~ Ky ) 


The integral 1^ in Eq. (10) is the Beta function. The integrals in Eq. (10) 
can be evaluated for known values of Y . 
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Gu was obtained similarly as, 


Gjj - lim [ RSA 2iY Ii + RSA' 2iY Ij. + RS I 2 + RS I 2 ] (11) 

A-+0 

Because R — -iP and S - -iQ, equation 11 reduces to 

G IX - lim [ -PQA 2iY Ii - PQA' 2iY Ii + PQ I 2 + PQ T 2 ] (12) 

A->0 

Equations (9) and (12) show that the mode I and mode II strain energy 
release rates depend on A . The terms A-^^ can be rewritten as 
e ±2i7 ln(A) As ^ approaches zero, the exponential functions have no well 
defined limits; they oscillate between +1 and -1. In the finite element 
analysis, this means that the computed mode I and mode II strain energy 
release rates will be dependent on the crack tip element size and do not 
show convergence as the crack tip elements are made smaller. 

The total strain energy release rate is the sum of Gj and Gjj (see Eq. 
7) and is 

G — Gj + Gjj 

- lim [ 2PQ I 2 + 2PQ I 2 ] 

A->0 

« [ 2PQ I 2 + 2PQI 2 ] (13) 

The terms involving A-^^ dct not appear in the total strain energy release 
rate. The total strain energy release rate is therefore independent of A. 
This analysis suggests that the total strain energy release rate will 
converge as the crack tip element size decreases but the individual mode 
components depend on A, and hence, do not converge. 
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On the other hand, if the two materials are so chosen that the 
oscillatory power, y , is equal to zero (i.e when + mb^A “ MB + A*a^a) » 

the mode I and mode II components in Eq. (9) and (12) do not depend on A 
and, thus, have a finite limit. Thus when Y-0, the individual as well as 
the total strain energy release rates must converge as the size of the crack 
tip elements decreases. 

Interfacial Crack in a Bi-material Orthotropic Plate 
Instead of the two isotropic plates considered earlier, consider a 
crack between two orthotropic materials. The singularity at the crack tip 
is also of the form -l/2±iY [7,11], The near field stress and displacement 
distributions have a similar form as Eqs . 4 and 6, but the expressions for 
P,Q,R, and S are too complex to write symbolically due to the complex 
functional forms. Therefore, the plane strain problem of two orthotropic 
materials subjected to remote tension was analyzed with a special crack tip 
singularity element which utilizes the near field solution that contains the 
oscillatory component at the crack tip [11]. From the numerical solution 
the constants P,Q,R, and S were evaluated for different combinations of 
orthotropic materials. The constants satisfied the following relations for 
all combinations: 

PQ + RS - 0 and PQ * RS (14) 

This means that like the isotropic case, for the bi-material orthotropic 
case the total strain energy release rates are independent of A (see Eqs (9) 
and (12)). Again, the individual components Gj and Gji will dependent on A, 
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because of the term in equations like 9 and 12. Thus the finite 

element solutions for the bi-material orthotropic case show trends like the 
bi -material isotropic case. 

The previous discussion is centered on plane problems. However, the 
focus problem is a Q3D problem. The Q3D problem can be thought of as a 
superposition of two problems. Problem 1 is with e x - € 0 and the rest of 
the five strains are zero. Problem 2 is with G x - 0 and the rest of the 
strains are non zero (see Eq. (1)). Problem 1 yields the non-homogeneous 
part of the solution, while problem 2 yields the homogeneous part. Thus one 
can concentrate on problem 2 to understand the behavior of the solution. 

But problem 2 is exactly equivalent to the plane strain problem. Thus the 
previously presented analysis for the plane problems is also valid for the 
Q3D problem under consideration. 


RESULTS AND DISCUSSION 

To evaluate the analysis presented in the previous section, an edge- 
delamination problem was studied using a laminate with two isotropic 
materials having zero or non-zero values of y . Next, the edge -delamination 
problem was analyzed using a resin-rich layer at the interface. The strain 
energy release rates for the resin layer case were compared with those from 
a 'bare' interface model to evaluate the accuracy of the 'bare' interface 
model . 
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Laminate with Isotropic Materials 

To demonstrate that the oscillatory part of the singularity is the 
cause of the non- convergence of the individual mode components, a laminate 
consisting of two isotropic materials A and B with ply thicknesses 3h and h, 
respectively, and a delamination at z - h was considered (i.e in Fig. 2(a), 
the top three plies are of material A and the bottom ply is of material B) . 
Materials A and B had the following properties: 


Material A : 

Ea 

= 10 X 10 6 

psi ; 

V A - 0.3 

Material B : 

Eb 

- 19.231 x 

10^ psi ; 

V B - 0.0 


These materials were chosen so that the y * 0 for plane strain conditions. 
Table 2 presents the individual as well as the total strain energy release 
rates when various size non- singular and quarter-point elements were used at 
the delamination tip. The finite element solutions yielded almost identical 
results for various element sizes at the delamination tip. 

Next, the material properties of the ply below the delamination was 
changed to 

Material B : Eg - 30 x 10 6 psi ; V B - 0.0 

From equation (5), this combination of materials gives a non-zero 
oscillatory power 0.087. Table 3 presents the strain energy release rates 
for this material combination with non- singular and quarter-point elements 
at the delamination tip. As expected, the individual modes do not 
convergence as the delamination tip elements are made smaller but the total 
strain energy release rate remains practically unchanged with mesh 
refinement. These results confirm that the oscillatory part of the 
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singularity is the cause of the non- convergence of the individual mode 
components . 


Laminates with Orthotropic materials 

As pointed out earlier, the stress singularity for delamination 
between two orthotropic materials is also of the form -1/2 ±i Y [7,11]. For 
the particular material properties considered and for the [0/135/90] s 
laminate with a delamination at the -35/90 interface, the singularity is 
-0.5 ± 0.0225i [7,16], (This value is determined by interpolation of values 
given in reference 7). The isotropic results presented earlier, indicates 
that the oscillatory part of the singularity is causing the non- convergence 
of the individual modes presented in Table 1. 

As pointed out in the introduction, in reality, a discrete or 'bare' 
interface does not exist between plies as was assumed earlier. When the 
laminates are cured, a thin resin layer develops between the neighboring 
plies. This resin layer is about 0.0004 in. thick. For convenience, the 
delamination was assumed to exist centrally within this resin layer. An 
eight ply laminate, [ 0/±35/r/90] s , with a delamination in the interface 
resin layer was analyzed. Note that the resin layer at interfaces other 
than the delaminated interface were not modeled because their influence on 
the strain energy release rates was considered negligibly small owing to 
their small thicknesses . As the delamination exists in a homogeneous 
isotropic material (the interface resin layer), the singularity has the 
classical square root power without the oscillatory component. Thus, the 
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individual components of the total strain energy release rate should 
converge as the size of the delamination tip elements are decreased. 

Three finite element idealizations with A/h r - 0.25, 0.125, and 
0.0625 were developed, where h r is the resin thickness (h r /h - 0.074). At 
the delamination tip, square non- singular elements were used. The finite 
element models had about 1800 nodes and about 5400 degrees of freedom. This 
large number of nodes was required because smooth transition was needed away 
from the delamination tip, and the delamination tip elements in the resin 
layer were much smaller than those used in the 'bare' interface models. 

Table 4 presents results obtained with these models. As expected, the 
individual as well as the total strain energy release rates showed very 
little change the for various values of A/h r values. Also, the Gm 
component is now positive. 

The 'bare' interface laminate was reanalyzed with the same physical 
delamination tip elements as the resin layer model. Table 5 presents the 
results with this model. Note that the model here has delamination tip 
elements that are an order of magnitude smaller than those presented in 
Table 1. These results show the same nonconvergent behavior as Table 1 
results. The results for the resin layer model (Table 4) and the 'bare' 
interface results of Tables 1 and 5 are plotted in Fig. 7. Note that A/h on 
the abscissa is plotted on the log scale. The horizontal dashed line 
represents the mode I and mode II components for the laminate with a thin 
resin layer at the delaminated interface. The results for the laminate with 
the 'bare' interface are not vastly different from those with the cocured 
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laminate if A/h values between 0.25 and 0.5 are used. However, larger 
discrepancies exist if much smaller A/h values are used. The maximum 
differences are less than 4 percent for mode I and about 2 percent for mode 
II if A/h values of 0.25 or 0.5 are used. Therefore, the finite element 
model with the 'bare' interface (A/h - 0.25 or 0.5) was a very good 
approximation to the case with the interface resin layer, although this 
model suffers from the non- convergence of the individual modes as the 
delamination tip elements are made smaller. This model is attractive 
because fewer elements are required and, hence, fewer degrees of freedom. 
The present results suggest that the size of the delamination tip elements 
should be one-quarter to one-half of the ply thickness. 


CONCLUSIONS 

The problem of edge -delamination of a long [0/±35/90] s composite 
laminate subjected to uniform axial strain was studied. This laminate had 
edge -delaminations at the -35/90 interfaces. Finite element models had 
either a ' bare' interface or a thin resin layer at the interface. A quasi 
three-dimensional finite element analysis was used to calculate the strain 
energy release rates. The virtual crack closure technique was used to 
obtain the individual components as well as the total strain energy release 
rates . 


The finite element analysis showed that the individual components 
of strain energy release rate did not converge when the ratio of the size of 
the delamination tip element to the ply thickness, (A/h), decreased for the 
'bare' interface case. In contrast, the total strain energy release rates 
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were unchanged as A/h decreased and agreed very well with the value 
calculated from a closed form equation derived from the rule of mixtures 
using classical laminated plate theory. This non- convergence was observed 
with non- singular , 8-noded parabolic elements and with collapsed quarter- 
point singularity elements at the delamination tip. However, the results 
obtained using the non- singular and singular elements agreed very well with 
each other. / 

Results for the resin layer case showed that the individual 
components as well as the total strain energy release rates remain 
practically unchanged as A/h decreased. Based on the studies performed, the 
following conclusions are drawn: 

1. The non- convergence of the individual components of the strain energy 
release rate calculated by quasi three-dimensional finite element analyses 
is due to the oscillatory part of the stress singularity. Continuum 
analyses show that the individual components will have no definite limit as 
the virtual crack closure size approaches zero. The total strain energy 
release rate, in contrast, has a well defined limit. 

2. When the materials were chosen such that the oscillatory component of 
the singularity is zero, the finite element solutions showed convergence for 
the individual components as well as for the total strain energy release 
rates as A/h was decreased This is in agreement with the continuum 
analyses . 
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3. When values of A/h were either 0.25 or 0.5, the individual components 
of the strain energy release rate obtained with the 'bare' interface model 
agreed well with those obtained with the resin layer model. The maximum 
difference between the two models was less than 4 percent for the mode I and 
2 percent for the mode II components. 


I 

i 
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APPENDIX 


RELATIONSHIP BETWEEN THE MATERIAL PROPERTIES 
FOR EQUAL CONTRACTIONS e y NEAR THE CRACK TIP 


The purpose of this appendix is obtain a relationship between the 
material properties of materials A and B such that the contractions in the 
y-direction very near the crack tip are identical (see Fig. 6). Consider 
the stresses a z and Oy in the two materials very near the crack tip and 
along the z — 0 line, 


a zk “ CT zB 
a yk “ a yB 


*1 



(Al) 


The corresponding strains are 

€ yA “ t a yA " v k a zk ' \ a xk ]/ E A (A2) 

and 

e yB “ l CT yB ' V B CT zB * V B ff xB 1/ E B (A3) 

where 

o x j - 0 for plane stress 

- v j ( Oyj + a z j ) for plane strain, (A4) 

and j— A,B. 


If the strains e y are same for the two materials A and B then by 
equating Eq. (A2) and (A3) one can show that 

MA + MB k A ” MB + MA k B (A5) 
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The and /ig are the shear moduli of materials A and B, respectively, and 
kj is given by, 

kj - 3 - 4 v j for plane strain and 

3 - vj (A6) 

- for plane stress 

1+V j 

and j- A,B. 

Thus if the material properties of the two materials are such that Eq. (A5) 
is satisfied then the stress-distribution very near the crack tip is given 
by Eq. (Al) and it does not contain the oscillatory part V . However, if 
the material properties do not satisfy Eq. (A5) then the simple square root 
distribution alone cannot maintain the compatibility along the bond line and 
very near the crack tip. 
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Table 1: Individual and total strain energy release rates for a 
[0/±35/90] s laminate : Non-singular and quarter-point 

singularity elements. 


Percentage of individual modes 

A/h - 

0.5 

0.25 

0.125 

0.0625 

Mode I 

21.22 

22.41 

(21.89)* 

23.49 

(23.38) 

24.65 

(24.69) 

Mode II 

79.98 

78.82 

(79.10) 

77.76 

(77.61) 

76.58 

(76.31) 

Mode III 

-1.20 

-1.23 

(-0.99) 

-1.25 

(-0.99) 

-1.24 

(-1.00) 

G 

0.0407 

0.0406 

(0.0402) 

0.0407 

(0.0402) 

0 . 0406 
(0.0401) 

e 2 o t En 


* Values in parentheses obtained with quarter-point singularity 
elements . 


[ G/(e 2 Q t E n )] CLT - 0.0407 
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Table 2: Individual and total strain energy release rates for a 

laminate with two isotropic materials: Non-singular and 

quarter-point singularity elements. 

E A - lOx 10*> psi ; V A “ 0.3 

E b - 19.231x 10 6 psi ; V B - 0.0 

7- 0.0 



Percentage 

of individual 

modes 

A/h - 

0.25 

0.125 

0.0625 

Mode I 

39.20 

39.21 

39.19 


(38.76)* 

(38.82) 

(38.78) 

Mode 11 

60.80 

60.70 

60.81 


(61.24) 

(61.18) 

(61.22) 

Mode III 

0.0 

0.0 

0.0 

G 


0.0137 

0.0137 

0.0137 

e 2 o t E A 

(0.0137) 

(0.0137) 

(0.0137) 


* Values in parentheses obtained with quarter-point singularity elements. 
[ G/(6 2 o t E a )] C lt - 0.0137 
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Table 3: Individual and total strain energy release rates for a 

laminate with two isotropic materials: Non-singular and quarter- 

point singularity elements. 


- 10 xlO** psi 


V A “ 0.3 


Eg - 30 xlO^ psi 


V B “ 0.0 


y - 0.087 



Percentage of individual 

modes for 

A/h - 

0.25 

0.125 

0.0625 

Mode I 

47.81 

49.44 

51.06 


(46.94)* 

(48.59) 

(50.19) 

Mode II 

52.19 

50.56 

49.94 


(53.06)* 

(51.41) 

(49.81) 

Mode III 

0.0 

0.0 

0.0 

G 

0.0177 

0.0177 

0.0177 

e2 o 

(0.0178) 

(0.0177) 

(0.0177) 


* Values in parentheses obtained with quarter-point singularity elements. 



Table 4: Individual and total strain energy release rates for a 
[0/±35/r/90] s laminate with delamination in the resin 
layer: Non- singular elements. 

(h - 0.0054 in.; h r - 0.0004 in.; h r /h - 0.074) 



Percentage of 

individual modes 


A/h r - 

0.25 

0.125 

0.0625 

Mode I 

18.65 

18.68 

18.68 

Mode II 

81.13 

81.10 

81.09 

Mode III 

0.22 

0.22 

0.23 

G 

S 2 o t e a 

0 . 0404 

0 . 0404 

0.0404 


[ G/ (e 2 0 t E a )] clt - 0.0404 ; Here t -8h. 


30 


Table 5: Individual and total strain energy release rates for a 
[0/±35/90] s laminate with delamination between -35 and 
90 plies: Non-singular elements. 



Percentage 

of individual modes 


A/h - 

0.0370 

0.0185 

0.00925 

Mode I 

25.53 

26.74 

28.00 

Mode II 

75.68 

74.43 

73.12 

Mode III 

-1.21 

-1.17 

-1.12 

G 

e 2 0 t E n 

0.0406 

0.0406 

0.0406 


[ G/ (e 2 0 t E n )] C LT * 0-0^07 
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Fig. 2: Typical x=constant plane and free body diagram 
showing equilibrium of top group of plies 





(a) Non-singular elements 
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Fig. 4: Nodes used in the strain energy release rate calculations 
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Fig.5: Comparison of strain energy release rates obtained with 
non-singular and quarter-point singularity elements 
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Fig. 6: Bi-material plate with 
an interfacial crack 


Results of the resin layer model 
'bare* interface model 
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Fig. 7: Comparison of strain energy release rates from 
the resin layer and 'bare* interface models 
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